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Abstract. Discovered over forty years ago, Gamma-Ray Bursts (GRBs) re- 
main a forefront topic in modern astrophysics. Perhaps the most fundamen- 
tal question associated with GRBs is the nature of the astrophysical agent (or 
agents) that ultimately powers them: the central engine. In this review, I fo- 
cus on the possible central engines of long-duration GRBs, and the constraints 
that present observations place on these models. Long GRBs are definitively 
associated with the deaths of massive stars, but whether the central engine is 
an accreting black hole or a rapidly-spinning, highly-magnetized neutron star 
(a "proto-magnetar" ) remains unsettled. This distinction has been brought into 
particular focus by recent MHD simulations of the core-collapse of massive, ro- 
tating "coUapsar progenitors," which suggest that powerful magneto-centrifugal 
outflows from the proto-neutron star may stave off black hole formation entirely. 
Although both black hole and magnetar GRB models remain viable, I argue that 
the magnetar model is more mature in the sense that it provides quantitative 
explanations for the durations, energies, Lorentz factors, and coUimation of long 
GRB outflows. Given these virtues, one promising strategy to break the present 
stalemate is to further develop the magnetar model until inescapable (and fal- 
sifiable) predictions emerge. This course of action signals a renewed challenge 
to translate time-dependent jet properties (power, magnetization, and Lorentz 
factor) into observables (gamma-ray light curves and spectra). 



1. Introduction 

Despite being discovered over forty years ago (Klebesadel et aL 1973), Gamma- 
Ray Bursts (GRBs) largely remain enigmatic. New observations, especially by 
the Swift and Fermi satellites, are continuing to advance our understanding of 
these events and many puzzles have been resolved. However, as new phenomenol- 
ogy has been uncovered, previous assumptions have been challenged and many 
new questions have emerged (see, e.g., Gehrels et al. 2009 for a recent review). 
Furthermore, answers to many of the most fundamental questions associated 
with GRBs remain (somewhat exasperatingly) elusive; for instance, "What emis- 
sion mechanism creates the gamma-rays? (e.g., synchrotron, inverse-Compton, 
hadronic?)," "What is the composition of the GRB-powering outflows? (e.g., 
baryons or e~/e^ pairs?)," "In what form is the outflow's energy stored (e.g., 
kinetic energy or Poynting flux?)," "How is this ordered energy 'randomized' and 
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how are the radiating charges accelerated?," "How are the magnetic fields that 
appear necessary for the prompt and afterglow emission produced? (e.g., are 
they generated by shocks or carried out from the outflow's source?)." 

Perhaps the key question associated with GRBs, and the focus of this re- 
view, is the nature of the astrophysical agent (or agents) that ultimately powers 
them: the "central engine" (CE). Soon after their discovery, there were per- 
haps more theories for GRBs than theorists (Ruderman 1975). The presently 
appreciated requirements of supernova-scale energies, short timescales (down to 
milliseconds), and relativistic speeds (Lorentz factors ^ 100) have, however, sig- 
nificantly narrowed the possibilities: GRBs are almost certainly the result of 
stellar-mass black holes (BHs) or neutron stars (NSs) being formed or undergo- 
ing catastrophic rearrangement (e.g., Katz 1997). All plausible CEs are powered 
by either accretion or rotation, which makes a deep gravitational potential and a 
reservoir of significant angular momentum the key ingredients of any model. A 
full understanding of GRBs may ultimately require comprehending the interplay 
between the physics of relativistic fluid dynamics, ultra-strong gravity, strong 
electromagnetic fields, nuclear /weak interactions, and plasma processes such as 
collisionless shock formation, non-thermal particle acceleration, and magnetic 
reconnection. The ways in which these nominally disparate physical processes 
conspire to produce a GRB makes studying the CE both exciting and uniquely 
challenging. 

Our discussion will focus on simple theoretical models of GRB CEs in order 
to elucidate both the abilities and the limitations of these prime movers. Fun- 
damental questions about GRBs (such as those raised in the first paragraph) 
are typically addressed by taking an agnostic approach towards the nature of 
the CE. Pursued in the spirit of "model independence," this approach also often 
allows the CE undue freedom, making it "a flexible source of power whose prop- 
erties are only limited by its total mass and the referees of theoretical papers" 
(Blandford 2002). It is important to emphasize that rather stringent constraints 
can be placed on GRB models by employing a self-consistent physical model 
for the CE. It is certainly true, for instance, that the formation, collimation, 
and stability of ultra-relativistic jets remains a formidable unsolved theoretical 
problem (e.g. McKinney &: Blandford 2009). However, other characteristics of 
the CE (such as the mass loading of the jet and energy budget) can, in some 
cases, be evaluated with more confidence. These aspects of the problem are 
more theoretically tractable largely because the immediate vicinity of the CE is 
often highly collisional and (local) kinetic equilibrium is assured; in this sense, 
the "extreme" environment of the CE can be a blessing. 

The review is organized as follows. I begin in (^ 12?]) with a brief historical 
background and summary of the basic observational constraints on CE models. 
It has long been known that GRBs come in (at least) two varieties ("long" 
and "short"), separated by their duration ~ 2 seconds (Mazets et al. 1981) 
and spectral "hardness" (Kouveliotou et al. 1993), and Swift has enabled the 
discovery that long and short GRBs indeed have distinct stellar progenitors 
(see, e.g., Berger 2009). In this review I focus on long-duration GRBs, which 
are generally thought to result from the deaths of massive stars (Woosley &: 
Bloom 2006). For a discussion of short GRB central engines, I refer the reader 
to excellent reviews by Nakar (2007) and Lee &: Ramirez-Ruiz (2007). The body 
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of the review divides into a discussion of collapsar (black hole; and proto- 
magnetar central engine models. I conclude and discuss the prospects for 
distinguishing these models in (^5?|). Since GRBs comprise an enormous research 
field, the topics covered here are necessarily limited and my reference list is only 
representative. Furthermore, the discussion and emphasis necessarily reflects 
my own work and biases. I apologize for any glaring oversights or omissions 
and refer the reader to several other excellent reviews for a different perspective: 
Lyutikov & Blandford (2003); Spruit (2004); Zhang & Meszaros (2004); Piran 
(2005); and Thompson (2008). 

2. Historical Development and Observational Constraints 

GRBs were discovered using the Vela satellites in the late 1960s, but they were 
not announced until 1973 (Klebesadel et al. 1973). GRBs manifest as a deluge of 
high-energy emissions with characteristic durations of ~ milliseconds to minutes; 
the light curves display rapid variability and a non-thermal (broken power-law) 
spectrum that peaks in the sub-MeV energy range but that often possesses a 
high-energy tail with significant power above ~ MeV (known as a "Band" spec- 
trum; see e.g. Fishman &: Meegan 1995 for a review). The fact that GRBs are 
variable on timescales typically <C than their total duration suggests that the 
outflow producing the GRB is best described as a quasi-continuous wind from 
the CE, rather than a singular explosion (Kobayashi et al 1997; see, however, 
Narayan & Kumar 2009). A recent discovery made with Fermi is that some 
bright GRBs are accompanied by an additional, typically flat-spectrum GeV 
component (Abdo et al. 2009b), which is generally delayed with respect to the 
normal Band-spectrum photons by a few seconds. While the ^ MeV flux from 
GRBs normally decreases rapidly once the GRB is finished (e.g. Barniol-Duran 
&; Kumar 2009), the GeV component appears to decay relatively gradually in 
time (oc t~^'^). Although the origin of this emission is currently debated, it 
may result from the onset of the afterglow emission as the GRB-producing out- 
flow first interacts with the surrounding, circumburst medium (e.g. Kumar &: 
Barniol-Duran 2009). 

As cosmological sources, the implied isotropic energies of GRBs are enor- 
mous, ranging from 10^^ ergs (for some short GRBs and local, sub-luminous 
long bursts) up to ^ 10^^ ergs ~ Mq(? in extreme cases (e.g. the recent Fermi- 
detected burst GRB 080916C had Ei^o ~ 8 x 10^^ ergs; Abdo et al. 2009a). If 
such a huge energy flux truly originated from a region with a size ~ 100 km 
(typical of that inferred from the observed millisecond variability), the opacity 
to pair-production in the high-energy photon power-law tail would be enormous 
(the implied optical depth for photon-photon and photon-electron interactions is 
r^_^,Te-'y ~ 10^^ for typical parameters). This requires the formation of a ther- 
mal pair photosphere, which is inconsistent with the observed non-thermal spec- 
trum. In order to overcome this "compactness problem," the GRB-producing 
region must be expanding towards the observer ultra-relativistically, with a bulk 
Lorentz factor F ^ 100 (e.g., Lithwick & Sari 2001). One reason that relativistic 
bulk motion alters the conclusion that high-energy gamma-rays may escape is 
that the observer, by naively inferring the size of the emitting region from the 
variability, is tricked by relativistic effects into believing that the emitting re- 
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gion is more compact than its true physical dimensions (indeed, current models 
place the location of GRB emission at radii ranging from 10^^ — 10^® cm, not 100 
km; compare, for instance, the emission models of Thompson 2006 and Lyutikov 
&; Blackman 2001). Relativistic speeds in GRB outflows have been confirmed 
observationally by the angular expansion rate inferred from the quenching of 
scintillation in a radio afterglow (Frail et al. 1997) and by direct VLBI imaging 
of an expanding GRB blast wave (Taylor et al. 2004). 

Since GRB outflows are ultra-relativistic, observers are only in causal con- 
tact with a small solid angle (~ 1 /r^) on the surface of the outflow (as subtended 
by the CE). Thus, even if GRB outflows are significantly collimatcd, evidence for 
jet-like angular structure only becomes apparent long after the main GRB event, 
once the flow has slowed to a Lorentz factor T ~ ^/dj, where 9j is the open- 
ing angle of the jet. Indeed, observational evidence such as "jet breaks" in the 
late-time afterglow (e.g. Rhoads 1997; Cenko et al. 2009) and the calorimetry 
provided by late-time radio afterglows (e.g. Berger et al. 2003) strongly suggest 
that GRB outflows are coUimated to some degree. This implies that although 
the isotropic-equivalcnt energies of GRBs are enormous, the true energy budget 
(after correcting for beaming) is much lower, typically comparable to the kinetic 
energy of a SN, ~ 10^^ ergs (Frail et al. 2001). What makes the GRB CE unique 
is therefore not necessarily its large energy budget, but rather its ability to place 
a significant fraction of this energy into material with ultra-relativistic speeds. 
This is an important diagnostic because it implies that the outflow remains rel- 
atively "clean" despite originating from what is potentially a rather dense and 
messy environment around the CE. A jet with energy ~ 10^^ ergs must, for 
instance, entrain ^ 1O~^M0 to achieve F ^ 100. Because all viable CE models 
are surrounded by signiflcantly more mass than this (~ 10~^ — IOMq), one of 
the major challenges of any model is to prodiice an outflow that avoids being 
polluted by too much of this surrounding material. 

An important discovery made with Swiffs X-ray telescope (XRT) is that 
GRBs are often followed by late-time X-ray flaring on timescales of ~ minutes- 
hours following the burst (e.g. Burrows et al. 2005). Because of the flares' rapid 
evolution and other similarities to prompt GRB emission, they most likely result 
from continuing CE activity (e.g. Lazatti & Perna 2007). The "plateau" phase 
observed in many X-ray afterglows on timescales ~ 10^ — 10^ seconds may also 
be powered by the CE (e.g. Panaitescu 2008; Kumar et al. 2008; see, however, 
Uhm &; Bcloborodov 2008). Thus, an important additional constraint on CE 
models is the ability to be active long after the GRB itself. 

3. Collapsar (Black Hole) Model 

The earliest cosmological CE models, such as neutron star mergers, focused 
(justiflably) on systems thought to naturally provide a relatively baryon-clean 
environment from which to launch a relativistic jet. In 1993 Stan Woosley made 
the bold suggestion that GRB jets may be produced even in the dense environ- 
ment of the core-collapse of a massive star (Woosley 1993). In this "collapsar" 
model, the GRB is powered by accretion of the stellar envelope onto a BH that 
forms soon after the collapse (e.g. Proga &: Begelman 2003). A jet produced 
by the accreting BH burrows through the collapsing star, producing a channel 
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through which the relativistic outflow can then escape (MacFadyen & Woosley 
1999; Zhang et al. 2004; Morsony et al. 2007). In order for a centrifugally- 
supported disk to form, the core of the stellar progenitor must itself be rapidly 
rotating; since only a small fraction of all massive stars are likely to satisfy 
this criterion (e.g., Cantiello et al. 2007; Langer et al. 2008), this provides an 
explanation for why GRBs are such a rare phenomena. The duration of the 
GRB in this model is presumably similar to the ~ 100 second in-fall time of the 
progenitor's He core. 

The original collapsar model envisioned a "failed" SN, in which most of the 
star accretes onto the BH instead of becoming unbound. Somewhat ironically, 
primary support for the general collapsar picture has come from the association 
of some long-duration GRBs with bright, "hyper-" energetic (E^n ~ 10^^ ergs) 
Type Ic SNe (e.g. Stanek et al. 2003; Modjaz et al. 2006), which originate from 
stellar progenitors that have lost their outer H/He envelopes (i.e., Wolf-Rayet 
stars). The association of long GRBs with regions of massive star formation 
(Bloom et al. 1999; Fruchter et al. 2006) has also been key to establishing that 
most (and possibly all) long-duration GRBs are produced by the core-collapse 
of massive stars (Woosley &: Bloom 2006). 

Since energetic SNe are in fact observed coincident with GRBs (i.e. they 
are hardly a complete "failure"), if the collapsar picture is correct this implies 
one or both of the following (cf. MacFadyen et al. 2001): 

1. The explosion mechanism associated with GRB-SNe is fundamentally dif- 
ferent than that associated with the death of "normal" (slower rotating) 
massive stars. MacFadyen &: Woosley (1999), for instance, propose that 
powerful winds from the accretion disk itself may explode the starQ It 
has been argued, however, that GRB-SNe are not markedly different that 
non-GRB Type Ic SNe in either their ^^Ni yield or peak brightness (e.g. 
Soderberg et al. 2006), although GRB-SNe are more energetic than "typ- 
ical" Type Ic SNe. 

2. The BH forms only after a delay, due to the "fallback" of material that re- 
mains gravitationally bound despite a successful SN (e.g. Chevalier 1993). 
In this case, the central compact object necessarily initially goes through 
a NS phase (which must have a decisive role in powering the SN), but 
the GRB is produced later, following BH formation. The delay until ac- 
cretion and collapse cannot be too long in this case, however, (i.e. within 
~ 10 — 100 seconds after core bounce) or the accretion rate onto the newly- 
formed BH would be too low to explain a typical GRB luminosity. 

Although long-duration GRBs are definitively associated with the deaths 
of massive stars, this does not establish that the CE is a BH. Single star evo- 
lutionary calculations with SN explosions put in by hand suggest that stars 
with initial masses greater than ~ 25Mq form BHs instead of NSs (Woosley 
&: Weaver 1995). However, this conclusion depends on the detailed SN mech- 
anism/energy and the stellar progenitor's pre-collapse structure, both of which 



^ A disk wind is probably not the mechanism for most core-collapse SNe, however, since accretion 
would spin-up the NS and the inferred rotation rates of pulsars at birth are typically much too 
low (e.g., Kaspi & Helfand 2002). 
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are poorly understood (see below). Indeed, our present observational under- 
standing of the mapping between high mass stars, the Wolf-Rayet progenitors 
of GRBs, and their compact-object progeny is far from complete (e.g., Smith & 
Owocki 2006). 

A satisfactory understanding of the origin of long-duration GRBs ultimately 
requires an understanding of the SN mechanism(s) as a function of progenitor 
mass, rotation, metallicity, and binarity. Modern core collapse simulations find 
that the outward-propagating shock produced at core bounce initially stalls at a 
radius ~ 200 km due to neutrino and photodissociation losses. It has long been 
suspected that neutrino heating from the proto-NS may be crucial to shock 
"revival" and a successful explosion (e.g. Bethe & Wilson 1985). However, 
spherically symmetric simulations fail to produce an explosion (e.g. Mezzacappa 
et al. 2001) and although more accurate and physical simulations are continually 
being implemented (e.g. Scheck et al. 2006) and multi-dimensional effects are 
likely a crucial ingredient (e.g. Murphy and Burrows 2009), it seems unlikely 
that neutrino heating alone is capable of powering the ~ lO^'^ erg hypernovae 
associated with long GRBs (e.g. Burrows et al. 2007). 

The massive stellar progenitors of GRBs are, however, far from typical. The 
requirements for the collapsar (and, indeed, any GRB model) are rapid rotation 
and a strong, large-scale magnetic field {'^ 10^^ G). The later is required to 
produce a stable relativistic jet (e.g. McKinney et al. 2006; Tchekovskoy et 
al. 2008)0 and rapid rotation and strong magnetic fields likely go hand-in-hand 
in core-collapse: differential rotation provides a source of free energy to power 
field growth, via e.g. an a — ^2 dynamo in the convective proto-NS (Duncan & 
Thompson 1992) or the magneto-rotational instability (MRI; e.g. Akiyama et 
al. 2003). 

The crucial question then arises: Do SNe indeed "fail" and lead to 
BH formation when their progenitor core is rapidly rotating (and, 
as a likely consequence, possesses a strong magnetic field following 
collapse)?" or stated more directly: Are the requisite initial conditions 
for the collapsar model self-consistent? 

The additional energy reservoir (rotation) and means for extracting it (mag- 
netic fields) make magneto-rotational effects a much more promising avenue for 
producing energetic "hypernovae" than neutrinos alone. While studies of mag- 
netic fields and rapid rotation in core collapse have a long and storied history 
(e.g. LeBlanc & Wilson 1970; Bisnovatyi-Kogan 1970; Symbalisty 1984; Wheeler 
et al. 2000), only recently have simulations captured both the effects of MHD and 
neutrino transport (e.g. Burrows et al. 2007; Dessart et al. 2008; hereafter DOS). 
D08 calculate the collapse of a rotating 350C progenitor model (35Mq ZAMS) 
and collapsar candidate of Woosley & Heger (2006) using a 2D axisymmetric 
MHD code with fiux-limited neutrino diffusion. They endow the pre-collapse 
core with a magnetic field that results in a ~ 10^^ G field strength when com- 
pressed to nuclear densities, thus mimicking the magnetic field expected from 
the saturated state of the MRI (their simulations do not resolve the MRI itself). 



^Neutrino annihilation along the rotation axis of the black hole is another possibility for powering 
a relativistic outflow (e.g. Zaldamea & Beloborodov 2008). However, the efficiency of this 
process is relatively limited, and it may have difficulty explaining very luminous GRBs. 
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Figure 1. Left Panel: Entropy map a.t t = 666 ms after core bounce, 
from the magnetized, rotating "Ml" core-collapse calculation of Dessart et 
al. (2008) using the 35Mq ZAMS collapsar progenitor model of Woosley & 
Heger (2006). Note the presence of a bipolar MHD outflow, simultaneous 
with continued accretion of the stellar envelope through a disk in the equa- 
tor. Right Panel: Total proto-NS mass Mns as a function of time after core 
bounce. A/ns initially grows due to accretion. However, for model Ml the 
(magneto-centrifugally-driven) mass loss at high latitudes eventually exceeds 
the accretion rate and Mjqg begins to decrease. BH formation is thus unlikely. 
Reproduced with the kind permission of L. Dessart and The Astrophysical 
Journal. 



A sample of D08's results are illustrated in Figured] (cf. Burrows et al. 2007). 
Soon after the collapse to nuclear densities, a bipolar MHD outflow develops from 
the newly- formed proto-NS. The explosion is not initially successful over all solid 
angles, as matter continues to accrete through a disk near the equator. However, 
the angular momentum accreted by the NS maintains its rapid rotation, which 
enhances the neutrino-driven mass-loss rate from mid-latitudes due to magneto- 
centrifugal "slinging" (e.g. Thompson, Chang, & Quataert 2004; Metzger et 
al. 2007). Importantly, in their strongly magnetized model (Ml), the mass-loss 
rate at high latitude eventually exceeds the accretion rate and for t ^ 300 ms 
the NS's mass begins decreasingl Although D08's simulations cannot address 
the possibility of later fall-back, and a different progenitor angular momentum 
profile could change the conclusions, their results are nonetheless suggestive: a 
core self- consistently endowed with the requisite properties to produce a GRB 
in the collapsar model (strong B and rapid rotation) may not result in a BH at 
allEI 

If a BH is not created following core collapse, then a rapidly spinning, 
highly magnetized proto-NS (a "proto-magnetar" ) likely remains behind in the 
cavity produced by the (successful) outgoing SN shock. This naturally leads 



■^Note that this conclusion does not imply that BHs cannot form from very massive stars. It 
does suggest that BHs may not form from the subset of rapidly-rotating progenitors which are 
responsible for hyper-energetic SNe and GRBs. 
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to a second possibility for producing a long-duration GRB: the "millisecond 
magnetar" model. 

4. Proto-Magnetar Model 

Usov (1992) proposed that GRBs may be powered by the spin-down of a ncwly- 
formcd, rapidly-spinning magnetar. The baseline motivations for the millisecond 
magnetar model are two-fold: (1) the rotational energy of a maximally-spinning 
NS (spin period P ^ I ms) of radius R is E^ot ~ 3 x 10^^{P/lms)-'^{R/Wkmy 
ergs, more than sufficient to explain the energetics of most GRBs; and (2) the 
vacuum/force- free electro-magnetic spin-down power of the NS for a assumed 
(equatorial) surface dipole field strength S^ip is approximately (Spitkovksy 2006) 

Eff ~ 5 — 1 ~ 10 ergss ; — , (1) 

similar to observed (beaming-corrected) GRB luminosities E ~ 10^^ — 10^^ ergs 
s-i for P ~ 1 ms and B^ip ~ 10^^ - 10^^ G (Blackman & Yi 1998; Wheeler et 
al. 2000). Though initially speculative because observed radio pulsars generally 
have Sdip ^ 10''^^ G (Manchester 2004) , evidence for the existence of magnetars 
has increased significantly in the past decade (Kouveliotou et al. 1998; Hurley 
et al. 2005). Indeed, Woods & Thompson (2006) argue that - 10% of Galactic 
NSs are formed with P^ip ^ 10^^ G, and the location of some magnetars within 
massive star clusters suggests that they may originate from massive stellar pro- 
genitors 40Mq; Muno et al. 2006). 

Magnetars are not, however, born in vacuum, nor are their magnetospheres 
necessarily well-described as force-free at early times. Like other NSs, magnetars 
are born hot and emit a substantial neutrino luminosity as they radiate their 
gravitational binding energy ~ 3 x 10^^ ergs (e.g. Burrows & Lattimer 1987). 
Despite its relatively weak interaction, this huge neutrino flux drives baryons 
from the surface of the magnetar (primarily via i^e/i^e absorption) during the first 
~ 10 — 100 seconds of its life (Thompson, Chang, & Quataert 2004; Metzger et 
al. 2007). Such a "neutrino- heated wind" is a generic consequence of a successful 
SN (e.g. Burrows et al. 1995). Indeed, Qian & Woosley (1996) show that the 
neutrino-driven mass-loss rate Mi, from a non-magnetized, non-rotating proto- 
NS is approximately given by 

M„ = ( (^r^ (^r, 

^ Vl052ergss-V V3MeV; \^1.4M0y VlOkmy ' 

where M, R, L^, and Ty are the proto-NS's mass, radius, electron neutrino 
luminosity, and neutrinosphere temperature, respectively. 

Metzger et al. (2007; hereafter M07) study the effects of magnetic fields and 
rotation on proto-NS winds by solving the equations of neutrino-heated MHD 
in the equatorial plane of the proto-NS (somewhat analogous to the Weber k, 
Davis 1967 model of the solar wind). By encorporating the cooling evolution 
of the proto-NS (e.g. L,j{t) and Ty{t) from Pons et al. 1999) and calibrating 
their results using multi-dimensional MHD simulations (Bucciantini et al. 2006) , 
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Figure 2. Energy loss rate E and magnetization {a = Fmax = E/Mc^) as a 
function of time since core bounce, calculated for a proto-magnetar with initial 
spin period P = 1 ms, fixed surface dipole field strength B^ip = 3 x 10^^ G, 
and three values of the NS mass Mns = 1-2Mq {dashed line), IAMq {dotted 
line), and 2.0MQ{solidline) (from the model of Metzger et al. 2007). This 
Figure illustrates that millisecond proto-magnetars produces outflows with 
the correct power {E ~ 10^° ergs s~^) and magnetization (maximum Lorentz 
factor Fmax ~ 10^ — 10"^) on the correct timescale (t ~ 10 — 100 seconds) to 
explain long-duration GRBs. The abrupt rise in cr at i ^ 30 — 100 second, 
which occurs when the proto-NS becomes optically-thin to neutrinos, provides 
a possible explanation for the observed sharp cut-off in GRB prompt emission 
(e.g. Barniol-Duran & Kumar 2009). 



M07 calculate how the mass-loss rate M and energy-loss rate E of magnetars 
evolve during the first ~ 100 seconds following their formation. From M{t) and 
E{t), they determine the "magnetization" a{t) = E/Mc^. For ultra-relativistic 
outflows (cj S> 1), CT represents the mctximiini potentictl Lorentz factor rmax? 

ob- 
tained if the outflow's Poyting flux is completely converted into kinetic power. 
Thus, a GRB-producing outflow must have a ^ 10^ — 10'^ to overcome compact- 
ness constraints { ^2^ - 

Figure [2] shows calculations of E{t) and a{t) from an updated version of 
M07's model for a proto-magnetar with B^ip = 3 x 10^^ and an initial spin period 
P = 1 ms, calculated for three different values of the NS mass (M = 1.2, 1.4, 
and 2.0 Mq). Initially E rises as the proto-NS contracts and spins-up, but E 
secularly decreases for t ;^ 2 s once the NS begins to spin down. At early times 

{t ~ few seconds) the proto-magnetar wind goes through a powerful {E ~ 10^^ 
ergs s~^) non-relativistic {a < 1) phase, which plays a crucial role in carving 
a bipolar channel through the stellar progenitor's envelope (see Fig. [3]), but is 
too mass- loaded to produce a GRB. At later times, however, a oc rises 

rapidly because M oc L^^^T^^^^ decreases quickly as the NS cools. The proto- 
magnetar wind becomes ultra-relativistic {a = Fmax ^ 100—1000) on a timescale 
t ~ 10 — 100 seconds, comparable to the duration of a typical GRB. 

We illustrate the robustness of this important result by estimating a explic- 
itly. Although at early times magneto-centrifugal enhancements to M^^ (eq. [2]) 
can be significant and Epp (eq. [1]) is not accurate for cr ^ 1, at later times 
(once the NS spins down a bit and a ^ 1), equations ([2]) and ([1]) are reasonable 
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approximations. The magnetization when a ^ 1 is therefore 



KW^^GJ V2ms/ UO^Oergss-V V2MeV/ \IAMq J \WkmJ 

where Rh = 271 Pc is the light cyhnder radius, the factor R^/2Ri^ is the fraction 
of the NS surface threaded by open magnetic flux for an aligned rotator (i.e. 
we assume only the "open zone" of the magnetosphere contributes to My), and 
we have scaled and Ty to typical values at t f» 30 seconds after core-bounce 
(Pons et al. 1999). Equation ([3]) shows that for the same magnetic fields and 
rotation rates necessary to power long GRBs (-Bdip ^ 10^^ G, P ~ 1 ms), the 
magnetar model naturally predicts the correct Lorentz factors (F ~ 10^ — 10^). 
Equation ([3]) is robust because My (eq. [3]) results from well-understood thermal 
physics and weak interactions at the proto-NS surface (Qian &: Woosley 1996; 
eq. [2]) and the basic electro-magnetic spin-down of a force- free rotator is now 
a solved theoretical problem (e.g. Spitkovksy 2006; eq. [1]). 

At later times (t ;^ 30 — 100 seconds) the proto-NS becomes optically thin 
to neutrinos and, consequentially, the neutrino-driven mass loss rate effectively 
vanishes. As this occurs a oc M~^ rapidly rises until finally plateauing at a 
very high value {a ^ 10^), once the mass- loss is dominated by electron/positron 
pairs driven by strong electric fields (e.g. Arons &: Scharlemann 1979) rather 
than neutrino heating. Since radiative processes likely become inefficient for 
ultra-high-o" flows (e.g. internal shocks move to very large radii and magnetic 
dissipation becomes slow; e.g. Giannios &: Spruit 2005), this z/-optically-thin 
transition plausibly signals the end of the prompt GRB phase in the magnetar 
model. The abrupt rise in cr at t 30 — 100 s may thus explain why prompt 
GRB emission is observed to cut off abruptly (e.g. Barniol-Duran &: Kumar 
2009). 

Although the above discussion illustrates that proto-magnetars produce 
outflows with the correct characteristics to explain long-duration GRB, it does 
not address how the outflow escapes from the surrounding, exploding star and 
becomes collimated into a bipolar jet. Using axisymmetric MHD calculations, 
Bucciantini et al. (2007, 2008, 2009) suggest a model for the coUimation of 
proto-magnetar outflows which is similar to that applied to explain the evo- 
lution and morphology of (much older) pulsar wind nebulae such as the Crab 
Nebula (e.g., Begelman &: Li 1992). The essential point is that the (exploding) 
star inertially confines the proto-magnetar wind (cf. Uzdensky &: MacFadyen 
2006; Komissarov & Barkov 2007), causing toroidal magnetic flux to accumu- 
late in a shocked nebula behind the outgoing SN shock. Magnetic hoop stresses 
confine the gas pressure in this "proto-magnetar nebula" along the polar axis. 
This in turn exerts an anisotropic pressure on the SN shell (the so-called "tube 
of toothpaste" effect), which rapidly opens a polar cavity that directs and col- 
limates the proto-magnetar wind into a bipolar jet. Snapshots of this process 
from Bucciantini et al. (2009) are shown in Figure [3l This represents the first 
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Case B: density (log) 




Case B: velocity 
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Figure 3. Evolution of the wind from a Bdip = 3xlO^^G,-P = l ms 
proto-magnetar, released inside an exploding 35Mq ZAMS star, from the 
axisymmetric MHD calculations of Bucciantini et al. (2009). The wind is 
injected at small radii with a time-dependent power E(t) and magnetization 
a{t) from the calculations of Metzger et al. (2007; see Fig. [2]). Left Panel: 
Density profile [log(g cm~'^)] at t — 2 seconds after core bounce. A magne- 
tized nebula develops behind the SN shock, whose anisotropic pressure begins 
carving a bipolar channel along the rotation axis. Right Panel: Velocity at 
t = 5 seconds in units of c (zoomed out). A complete channel opens along the 
polar axis. The nebular magnetic field redirects the proto-magnetar's outflow 
into an ultra-relativistic, bipolar jet with an opening angle 0jct 5 — 10°. 



self-consistent GRB jet calculation that extends from the central engine to the 
stellar surface and which is evolved for a time At ~ 10 s comparable to the 
duration of a GRB. Importantly, these calculations illustrate that once the jet 
breaks out of the star, its mass and energy flux map nearly one-to-one onto 
those injected by the proto-magnetar wind at small radii. This justifies using 
E{t) and a{t) calculated using a free wind model (such as those in Fig. [2]) as 
direct input into a GRB emission model. 



5. Conclusions and Future Directions 

Observationally distinguishing between BH and magnetar models for long-duration 
GRBs is a difficult task. A comparison between the models is summarized in 
Table [TJ In principle, a detailed comparison could be performed between theo- 
retical models and the observed prompt emission. Figure U] shows a prediction 
for the light curve expected from "isolated" magnetar birth (i.e. neglecting the 
effects of an overlying stellar envelope), assuming that radiation is generated 
as the free energy in the outflow is dissipated through internal shocks. The 
qualitative shape of the light curve resembles the time-averaged "envelope" of a 
surprisingly large fraction of long GRBs. Unfortunately, making definitive pre- 
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Table 1. Black Hole vs. Proto-Magnetar Long-Duration GRB Models 
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Figure 4. Luminosity of internal shock emission in proto-magnetar winds 
for a NS with initial rotation period P — 1 ms and three surface dipole 
magnetic field strengths: Bdip — 10^^ G {dashed line), 3x lO^'^ G {dotted line), 
and 10^^ G {solid line). The gradual onset of the emission begins once the 
magnetization of the outflow a reaches ^10 due to the large Thomson optical 
depth, which decreases as the outflow expands. The light curve was calculated 
for "isolated" magnetar birth (e.g. a quasi- isotropic outflow). Correcting for 
the collimation due to the surrounding stellar envelope (Fig. [S]), the observed 
luminosity is larger by a factor of ^ ~ 10^. The late-time BAT light 
curve from GRB 060614 (Gehrels et al. 2006), shown with a light solid line, is 
reproduced in a time-averaged sense by the Bdip = 3 x 10^^ G model. From 
Metzger, Quataert, & Thompson (2008). 
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dictions, even when the initial properties of the jet are known, is hindered by 
our ignorance of the dissipation and radiation mechanisms in the outflow (recall 
the unanswered questions from §1). 

The magnetar model does, however, make a few robust (albeit indirect) 
predictions. For instance, (1) at t ~ 10 seconds after core-collapse, outflows 
from higher luminosity bursts should possess higher Lorentz factors: 
although the spin-down luminosity is larger for more rapidly spinning, highly 
magnetized NSs, the mass-loss rate at this epoch is probably similar for all 
NSs (see eqs. [2] — [3]). Furthermore, (2) the Lorentz factor should be 
a monotonically-increasing function of time during the burst (Fig. [2]). 
The magnetar model also predicts that (3) for the first 10 — 100 s after core 
bounce (when the neutrino luminosity is large), the GRB-producing 
outflow should be baryon-dominated, unlike the pair-dominated outflows 
expected from radio pulsars (Arons &: Scharlemann 1979) or along flow-lines that 
thread the event horizon of a BH. If a baryon-rich outflow can be definitively 
ruled out on this early timescale (e.g. via searches for neutrino emission due to 
hadronic processes with, e.g., IceCube; Abbasi et al.2009), the magnetar model 
could be refuted. 

Unfortunately, similar predictions are not yet available for BH models, in 
part because it is still not clear how the mass and energy fluxes in jets from 
accretion disks are determined. Indeed, the jet's mass-loading may depend on 
uncertain diffusive processes from the jet walls (e.g. Levinson & Eichler 2000; 
McKinney 2005). Simply put, more robust predictions are possible from the 
proto-magnetar model because spheres are simpler than disks and the surface of 
a proto-NS is in kinetic equilibrium. Given these virtues, one promising strategy 
to break the present BH-magnetar stalemate is to transform the predictions for 
E{t) and a{t) from the magnetar model (Fig. [2]) into observable gamma-ray 
light curves and spectra (e.g. Fig. d]). Hopefully, then, robust and falsifiable 
predictions will emerge, and we may be able to discern the central engine by the 
process of elimination. This course of action signals a renewed theoretical and 
observational challenge to determine the source and mechanism of prompt GRB 
emission. 
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